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^ ' Abstract 

D ! 

The direct constraint on the parameters of short range pseudomagnetic 
interaction of free neutron with matter is obtained from the recent test 
' experiment on a search for neutron EDM by crystal-diffraction method 

(~| [I]. It is shown that this constraint on a product of scalar to pseudo-scalar 

coupling constants QaQp is better than that of any other method for the 
range A < 10~^cm. 
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^ ■ 1 Introduction 

■ Over last years a possibility to look for new hypothetic particles which results 
in a new short range Yukawa-type potential of fermion-fermion interaction is 
actively discussed. The spin-dependent short-range interactions may be induced 

OO ', by light, pseudoscalar bosons such as the axion invented to solve the strong CP 

■ problem This interaction is usually parameterised as [2] 



where gs and gp are nondimensional parameters of the scalar and pseudo-scalar 
coupling constants between the neutron and exchanged boson. A = h/niAC 
is typical parameter of the range of forces (Compton wavelength of axion). 
There are proposals to search this new type of interaction using gravitationally 
bound quantum states of a free neutron [3] and using a spin precession of the 
trapped ultracold neutrons in vicinity of bulk matter ^ . Both of these methods 
have a suitable sensitivity for the range of A > lO^^cm, but their sensitivity is 
extremely decreased for the range A < 10~^cm. 

Here we consider a possibility to use a neutron diffraction in the perfect 
non-centrosymmetric crystal to search a new type of short range interaction for 
the 10-1° < A < lO-^cm. 
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2 Diffraction in a non-centrosymmetric crystal 



Neutron diffraction in a non-centrosymmetric crystal was widely discussed within 
the framework of the project to search for the neutron electric dipole moment 
(nEDM) by the diffraction method 

Any crystal potential (nuclear, electric, new short range potential, ..) can 
be presented as sum of the potentials of different atoms placed into the crystal 
cell. For the periodic crystal structure it is convenient to present such potential 
as Fourier series over the reciprocal lattice vectors g 



V{r) = Vair -ra)^Y. ^9^'^"" = '^o + E 2i;g cos(gr + 0,,), (2) 
a g g 

where Va{r — ra) is the potential of single atom, Va is the atom position, Vg = 
Vg eKp(i(f)g), g = 2TT/d, d is the interplanar distance. Here we take into account 
Vg — Vlg, because we consider the real value potentials. 

g-harmonics of potentials can be found from the equation 



Vg 



In the case of nuclear potential 

27rft2 
mVr 



J d^r e~'9rv{r), (3) 



Vg = -—FT^S' (4) 



here m is the neutron mass, Vc is volume of crystal unit cell, Fg is the structure 
amplitude 

Fg^J2^~'^"'M9)e-'3^^. (5) 

i 

Here we sum over the atoms of unit cell, f{g) is the scattering amplitude of i 
atom. Wig is the Debye- Waller temperature factor. 

For the case of a non-centrosymmetric crystal the different potentials can 
be shifted to each others, by the other words, the phases (j)g of different crystal 
potentials can be not equal. For the case of electric potential this shift results in 
a large electric field affected the neutron in non-centrosymmetric crystal [71 [5] . 
We take the phase of the nuclear potential equal to zero so the value of the 
electric field affected the neutron will be 

Eir) = -grad (r) = 2vfg sin(gr + <pf). (6) 

where and 0^ are amplitude and phase of the g-harmonics of crystal electric 
potential accordingly. Let's consider the monopole-dipole interaction ([1]). Direct 
calculation of g- harmonic of Vsp{r) from ([3]) gives 



" " 2m 14 1 + g 
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where rig = g/g, F^^ and '^^^ are the amphtude and phase of structure factor 
fg^ of the crystallographic plane, fg^ is determined by the following sum 

i 

here Ai and is the mass and position of a corresponding atom in elementary 
cell. 

Neutron wave function is determined by the nuclear interaction, and for the 
case of close to Bragg direction passage through the crystal can be written as 



^(r) = e''^^ + —2 e^^f^ = e'^^ 

Ek — Ek„ 



jjN 

2A/ 



(9) 



where we take into account that the phase of nuclear harmonics equal to 
zero, kg^k + g, Ek = Tt^e/im, Ek^ = h^kl/2m, = h^U^ /2m, Ag = 
{kg — k'^)/2 is the parameter of deviation from the exact Bragg condition. 

Therefore, Vsp potential affecting the neutron in the crystal will be 

jjN 

Vsp = {i^{v)\Vsp{vmv)) = ^|F/^|sin$f = 



Ag 9 2mVc l + .g2A2 



We should note that for centrosymmetric crystal $g ^ = and in this case 
the mean potential affecting the neutron will be zero. One can see also from 
pUj) that this '"pseudomagnetic" ' potential is proportional to the parameter 
As = Ug/Ag determined by the deviation from the Bragg condition. That 
allows to control the value and sign of the potential selecting the neutrons with 
slightly different energies from the Bragg one. 

Interaction with such a potential will lead to the neutron spin rotation 
around the reciprocal lattice vector g by the angle 

V>sp = — (11) 
where r is the time of neutron travel through the crystal. 



3 The method sensitivity 

For example let's consider (110) plane of non-centrosymmetric quartz crystal 
and As = 0.5. For (110) plane g = 2.56 • lO^cm"!, F^^ = 51, sin($f ^) = 0.41, 
Vc = 113A3. The angle of spin rotation due to considered potential will be 

^5P = 0.36.1024[cm-3].^ff-^L (12) 
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where L is the crystal length. For the cold neutron beam at high flux reactor 
the measurement accuracy (j{(fsp) ^ 2 • 10~^ can be reached for 100 day of 
the statistic accumulation [I]. That allows to give the constraint for monopole- 
dipole interaction 

5S5P < lO-'Mcm'] • (ff' + l/A^) (13) 

for the L=50 cm. 

Recently the test experiment for the search for neutron EDM by crystal- 
diffraction method was carried out [Ij. This result already allows to give the 
direct constraint on a value of gs9p better than any other method for the A < 
10~^cm, see FiglH curve (1). 

The comparison of different constraints on gsgp is shown in Fig. [T] 
One should note that both neutron EDM interaction with crystal electric 
field and the spin-dependent short-range interaction lead to a neutron spin ro- 
tation about reciprocal lattice vector g. Therefore, the considered shot-range 
interaction can give a false effect for the neutron EDM experiment and vice 
versa. However, these two interactions will be different for different crystal- 
lographic planes, so in the case of nonzero effect they can be separated using 
different planes for measurement. 




X, cm 



Figure 1: Constraints on a value of coupling constants product gsgp- Curve (1) is 
the constraint from the crystal-diffraction nEDM experiment [1] (this work) and 
(2) is possible improvement of this method, (3) is gravitational level experiment 
[3], (4) is the UCN depolarization [10], (5) is proposal [4], (6) and (7) are the 
predictions of axion model with ~ 1 and 9 ^ 10~^° correspondingly [IO),[TT] 
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4 Conclusion 



Direct constraint on amplitude of T-odd monopole-dipole interaction of neutron 
with the matter was obtained. It is shown that the product of scalar to pseudo- 
scalar coupling constant gsgp < 10^^^ for the range 10^* < A < 10~^cm. This 
value can be improved on about 10'^ times for the full scale setup for the neutron 
EDM search by crystal-diffraction method, which is under construction now. 
This work is supported by grant RFBR-09-02-00446. 
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